Optical probing of shocks driven into overdense plasmas by laser hole-boring 
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Observations of the interaction of an intense Ao ~ 10 /im laser pulse with near-critical overdense 
plasmas (n e = 1.8- 3 n c ) are presented. For the first time, transverse optical probing is used to 
show a recession of the front surface caused by radiation pressure driven hole-boring by the laser 
pulse with an initial velocity > 10 6 ms -1 , and the resulting collisionless shocks. The collisionless 
shock propagates through the plasma, dissipates into an ion-acoustic solitary wave, and eventually 
becomes collisional as it slows further. These conclusions are supported by PIC simulations which 
show that the initial evolution is dominated by collisionless mechanisms. 



Collisionless shocks are associated with numerous nat- 
ural phenomena from the bow-shock ahead of the earth's 
magnetosphere [l] to the blast wave in supernova explo- 
sions [2]. They can accelerate particles to high energies, 
and thus are implicated in the generation of cosmic rays 
Q . It has been suggested that collisionless shocks can be 
similarly used to create energetic particle beams in the 
laboratory. In particular, the radiation pressure of pow- 
erful lasers has been proposed as a shock-initiating piston 
0-0 • A number of measurements of energetic ion beams 
have been attributed to acceleration by the space charge 
fields created at the target front surface when an intense 
laser strikes a target with initial density greatly exceed- 
ing the critical density, n ^> n c = eom e cc>o/e 2 , where ujo 
is the laser angular frequency [7[. 

In recent experiments, monoenergetic ion beams have 
been reported from targets of much lower (gaseous) den- 
sity. This was achieved by driving the target with longer 
wavelength (Ao ~ 10 fim) CO2 lasers jH, 0|. Proton 
beams with energies from ~ 1 to 20 MeV were observed 
with energy spreads < 5%. The acceleration was asso- 
ciated with shock-like structures moving away from the 
front surface of the plasma driven by the interaction with 
the laser. Furthermore, in 18] the most energetic beams 
were observed for n approaching n c . This is expected 
as momentum conservation implies the target surface's 
recession velocity should scale with mass density p as 
Vhb — v(l + v)I/P c i where 77 is the target reflecti vity 
due to the radiation pressure of a laser of intensity J[10|. 

The target front surface recession is called hole-boring 
and leads to a shock being driven into the target. 
Hole-boring has been experimentally identified from the 
Dop pler shift of light back reflected from the interaction 
ll| and its harmonics 



12|, and observed using x-ray 



refractometry [13[. There have also been a number of 
studies of collisionless shocks generated by intense lasers, 
either directly throu gh p onderomotive expulsion Q or 
from ablated plasma [15| , but these studies are generally 
limited to the under dense regime. Both x-ray studies and 



proton probing offer the ability to study to higher density 
but cannot directly measure density inside solid targets, 
where scattering is strong. 

This letter details a study of the shock driven into a 
just-over dense plasma by laser hole-boring. The plasma 
is driven by an infrared (CO2) laser, for which n c « 
10 19 cm -3 ; > 100 times lower than that for optical wave- 
lengths. This allows direct optical probing of the plasma. 
By varying the timing of the probe relative to the drive 
pulse, the shock evolution was studied, and three sepa- 
rate phases were identified. Initially the laser exerts a 
pressure on the opaque plasma, causing hole-boring and 
formation of either an electrostatic double layer at the 
surface, or a collisionless shock moving into the plasma, 
accelerating background ions. Once the laser stops, parti- 
cle reflection from the propagating shock structure causes 
it to decelerate rapidly, transitioning into a solitary ion- 
acoustic wave which continues to slow in the cooling 
plasma. As the shock speed decreases, the ion mean free 
path also decreases and the shock eventually becomes col- 
lisional. These findings are supported by particle-in-cell 
(PIC) simulations, which also help elucidate the mecha- 
nism of narrow energy spread ion beam generation. 

The experiment was performed with the Ao ~ 10 jam 
CO2 laser at the Accelerator Test Facility at Brookhaven 
National Laboratory. The laser was circularly polarised. 
Streak images of the laser time profile showed a train of 
pulses, caused by spectral modulation in the amplifiers, 
each of pulse-length r ~ 6ps (FWHM) and 18 ps pulse 
separation. On average 70% of the energy was in the 
leading two pulses, with a total integrated energy of 2- 
3 J. The pulse was focused with an f/3 off-axis parabolic 
mirror to a spot size wq w 70 //m (FWHM), giving a peak 
intensity I^5x 10 15 Wcm -2 and normalised vector po- 
tential ao ~ 0.5. The laser was focused ~ 0.7mm above 
a 2 mm diameter nozzle gas jet producing a helium gas 
neutral density profile, characterised by transverse inter- 
ferometry [16], approximated by an isosceles trapezoid of 
plateau width 1200 fim and side length 600 fim. 
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FIG. 1. a) Typical shadowgraphy and b) interferometry im- 
ages for t = 30 ps, n = 3n c . Vertical yellow lines indicate the 
initial critical density position. Laser enters from the left. 



The probe, a frequency-doubled (532 nm) Nd:YAG 
laser beam of r ^ 6 ps, was absolutely synchronised with 
the CO2 pulse. The timing between them was varied us- 
ing an optical delay line. The probe passed orthogonally 
through the plasma before being split; with one part im- 
aged directly for shadowgraphy, and the other going to 
a Mach-Zender interferometer, to diagnose plasma den- 
sity. Examples are shown in Fig. [1] for a peak initial 
density of n e = 3n c , at t = 30 ps after the start of the 
interaction, which is just after the second major pulse in 
the pulse train has passed. Bright features in the shad- 
owgram (Fig. QJi) highlight regions of rapidly changing 
refractive index and therefore density. The shadowgram 
shows a feature at the target front surface of ~ 100 /im 
transverse size narrowing to ~ 50 jam. Interferometry re- 
veals that this is a high density feature with a region 
of considerably lower density inside it (fig. lb), which 
is the channel formed by the initial 75 /im laser focal 
spot boring into the plasma. The feature ends ~ 45 jam 
longitudinally from the estimated initial critical surface, 
calculated from the gas density profile and supported 
by the off-axis plasma density profile. Assuming the 
front originated at the critical surface implies an average 
v = 1.5 x 10 6 ms _1 (±30%), which agrees reasonably with 
Vhb — 1-9 x 10 6 ms _1 , obtained using the laser parame- 
ters for this shot and taking p = mn e n c . This calculation 
ignores self focusing, the pulse train, and thermal effects. 

Fig. [2] shows a series of shadowgrams taken with vary- 
ing probe timings up to t = 1600 ps. The front seen in 
Fig. [Ik remains clearly defined, moving further into the 
plasma as time progresses. Horizontal features can be 
seen further into the plasma, particularly for t = 200 ps 
(fig. [2^). These are likely to have been created as the 
target was being ionised, since the gas target is initially 
transparent to the ionising laser. Below the shadowgra- 
phy images are the corresponding density maps derived 
from interferometry (fig. [2]i-f). The phase shift was as- 
sumed to be symmetrical around the laser axis, allowing 
Abel inversion to retrieve the density profile. The accu- 
racy of the recovered density maps is limited by asymme- 
try in the inter ferograms, caused by, for example, slight 



vertical density gradients and laser non-uniformities. 

Once sufficiently ionised, the plasma becomes opaque 
(n e > n c ) and the laser is reflected at the critical surface. 
It then exerts radiation pressure on the plasma, launching 
a shock into the initially relatively cold plasma. The peak 
density associated with this moving front exceeds the ini- 
tial background plasma density (fig. [2]b,c), and shows a 
steepening in the feature profile even at late times; clear 
indication of shock formation. Two probe timing scans 
were taken for peak densities n e G (1.8, 3) n c . The dis- 
tance of shock propagation from the initial critical sur- 
face position thus obtained is plotted in fig. [3^ as a func- 
tion of time. The shock initially moves rapidly into the 
plasma before slowing at later times, with the shock ve- 
locity initially higher for lower density. The mean veloc- 
ities for different parts of the shock trajectory are indi- 
cated on the plot. The previously estimated hole-boring 
velocity (vhb ~ 1-5 x 10 6 ms _1 ) is indicated, slowing to 
v sh « 5xl0 4 ms- 1 fort = 500- 1700 ps, for both n e ~ 1.8 
and 3n c . The ion mean- free-path due to (ion- ion) colli- 
sions for ions with v = Vhb is £mf ~ 30 mm [17[. This 
is much larger than the scale of the measured features 
(< 100 /am) and even the gas jet, indicating that the 
shock must initially be collisionless. At the latest times 
measured, as the shock slows, £ m /(oc v 4 ) decreases sig- 
nificantly, and collisions play a role in dissipation. 

To further investigate this behaviour, simulations 
of the interaction were performed using particle-in-cell 
(PIC) codes. 2D3V simulations were performed with the 
code OSIRIS [18]. A circularly polarised A = 10.6 jam 
pulse, with r = 6 ps, ao = 0.55, wl = 70 jam was focused 
on a fully ionised, initially cold (T e = Ti = 0) He plasma 
of density n e = 1.8- 6 n c with a step density gradient. 

Fig. 2^ shows the plasma density profile from the sim- 
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FIG. 2. (colour online) Time series of shadowgraphy (top) 
with corresponding density maps (below) for a,d) t = 200, b,e) 
500 and c,f) 1600 ps. The density is in units n c . Laser enters 
from the left. The shadowgraphy images have been Fourier 
filtered for clarity. The initial critical surface is indicated with 
a vertical line for all images. For all shots n e — 1.8 n c - 
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FIG. 3. (colour online) a) Shock propagation distance z s h 
with time from shadowgrams for n e £ (1.8, 3) n c (red circles, 
blue squares). Green solid line indicates initial Vhb, and black 
(dashed, dotted) lines represent the mean final v s h for n e £ 
(1.8, 3) n c ; b) from 2DPIC for n £ (1.8, 3, 6) n c (red circles, 
blue squares, green diamonds). The laser stops at t w 7.5 ps, 
indicated by vertical dotted line 



ulation just after the laser pulse. It shows a density in- 
crease at the front along with the cavity formed by the 
hole-boring, in agreement with what is seen experimen- 
tally. The front is formed by radiation pressure displac- 
ing plasma electrons to a collisionless skin-depth. The 
ions, which due to their greater mass are unaffected by 
the laser fields, accumulate at the front due to the space- 
charge field of the ensuing electrostatic double layer [l9[ . 
The large electrostatic field associated with this front also 
accelerates upstream ions. This leads to the formation of 
a large population of energetic ions that can be seen to 
outrun the shock front in fig. [4^. 

Since the plasma ahead of the front is initially cold, 
hole-boring causes a compression of the plasma. The 
piston-like effect of the radiation pressure launches a 
shock into the plasma. The motion of the shock front 
during the simulation is shown in fig. [3b for three differ- 
ent starting n e . Although initialised cold, the plasma is 
heated by vacuum heating [20|, due to the buckling of 
the front surface from hole-boring. Hence for lower n e , 
where hole-boring is faster, the rate of heating increases 
more rapidly, and the shock moves faster. 

To obtain high resolution, the simulation box for 2D 
simulations was limited to 60 /im longitudinally. As a re- 
sult a sheath quickly forms at the rear of the simulation 
box due to high energy electrons that leave the box. This 
sheath reflects electrons back into the plasma much more 
efficiently than in the experiment, where the plasma size 
was > 500 /im. This results in an unphysically high back- 
ground temperature T e , and consequently poor modelling 
of the shock front behaviour. To model the shock evo- 
lution with a more realistic box size, a series of ID PIC 
simulations were performed with the code EPOCH Id. A 
fully ionised He plasma with a top hat profile of n e = 2n c 
was irradiated with an ao = 1 circularly polarised laser, 
tl = 5 ps (gaussian) pulse with A = 10.6 /im. The simu- 
lation box was 200 /im, with 5000 cells and 2.5 x 10 5 par- 
ticles per species. The ions were set to be initially cold 



(Ti = eV). ID simulations cannot model the buckling 
of the plasma due to hole-boring. Therefore, to assess the 
influence of heating, the initial electron temperature was 
varied between T ei = 0-50 keV. These simulations are 
comparable to those in [4], as the lower intensity here is 
compensated by reduced densities to give similar piston 
velocities, and also in which discusses the role of ion 
temperature on the shock structure. 

Whilst the laser is on, the laser bores a hole into the 
plasma irrespective of T e $, although the hole-boring speed 
decreases slightly for higher temperatures due to the in- 
crease in the electron thermal pressure. For the highest 
T ei (= 50 keV), where Vhb is only just higher than the 
initial sound speed, c s = ^/kT e /mi, a series of solitary 
waves are pushed into the plasma. Decreasing the initial 
temperature of the simulation to 5 keV increases the ra- 
tio of piston velocity to the local sound speed, and the 
ion density profiles (e.g. fig. Hfc) show the formation of 
a collisionless electrostatic shock. The shock moves at 
v s h ~4x 10 6 ms _1 , compared with Vhb ~ 3 x 10 6 ms _1 , 
and so moves ahead of the critical surface. The density 
at this shock can be as high as n e = 10 n c . Fig. |4Jd shows 
that the shock width (~ c/uj p ) is less than 1 /im, which is 
why though it is clearly visible from shadowgraphy (0 a - 
c), it is not properly resolved in the interferograms ((2]d- 
f). The potential of this front (0.15 MV) is so large that 
the associated electric field, which reaches > 200 GVm -1 , 
efficiently reflects nearly all upstream ions. Phase space 
diagrams (fig. Hh) show that the reflected ions gain a 




FIG. 4. (colour online) a) Density distribution from 2DPIC 
for rte — 1.8n C r, Tei = at t = 9 ps. b) Evolution of ID 
He 2+ charge density for n e — 2n cr , T e i — 5keV, at t = 5 ps 
(blue, corresponds to end of laser pulse), 15 ps (red) and 30 
ps (green). v x x phase space for He 2+ from 1DPIC at the peak 
of laser for c) T e i = 5 and d) T e i = 1 keV. 



speed of 2v s h by 'bouncing-ofP the shock front. Whilst 
the laser is active, it is v s h that determines the acceler- 
ated ion velocity (not Vhb)- This results in a remarkably 
narrow energy spread for the reflected ions With 
T e i = 5 keV, there is partial reflection and partial trans- 
mission through the shock front. The transmitted ions 
display the characteristic for a shock solution with par- 
tial particle reflection, displaying numerous trailing per- 
turbations in the potential, as is evident in Fig. Hh (and 
UK) [21]. Further reducing T ei < 1 keV (so v hb > c a , 
fig. [4]i), a region of high ion density is formed right at 
the hole-boring front, creating a double layer moving at 
v s h ~ Vhb, which reflects nearly 100% of upstream ions. 

Therefore, the simulations show three separate inter- 
action regimes for T e ^> Tf. 1) for Vhb ^ 3c s , the laser 
launches solitary waves into the plasma from which there 
is no particle reflection; 2) for ^ 3c s , a shock forms 
moving ahead of the hole-boring front with partial parti- 
cle reflection; and 3) for Vhb ^> c s there is 100% reflection 
and v s h ~ Vhb 0, HH . When considering the interaction 
as a laser plasma ion source, there is therefore a tradeoff 
between higher maximum energies from a hotter plasma, 
where a small fraction of the particles are accelerated 
up to higher energies, and higher efficiency from a colder 
plasma where all the particles are accelerated to a smaller 
maximum energy. This explains the larger flux of accel- 
erated protons observed with a short, circularly-polarised 
laser beam [8], as opposed to a longer, linearly-polarised 
driver where heating is stronger 

At the end of the laser pulse, the potential associated 
with the shock is still sufficient to reflect ions. How- 
ever, as there is no longer a pressure being exerted by 
the laser, this particle reflection reduces the energy of 
the ion front, causing it to slow until reaching a criti- 
cal Mach number, M s , where particle reflection practi- 
cally stops. Assuming T{ = 0, one expects a lower limit 
of M s = 1.6 given an isothermal approximation for the 
plasma, whereas assuming maximum electron t rapp ing 



behind the shock gives an upper value of M s ~ 3 [2l| . In 
the simulations, at the end of the laser pulse the shock 
speed drops to v s h °° 3M s c s . For the 'cold' simulations 
(T e i = keV), this results in the shock practically stop- 
ping, not matching our observations. 

For simulations with a T e i ^ 0, T e decreases with time 
as energetic electrons leave the simulation box, causing 
the soliton speed to exceed M s c s . This stimulates fur- 
ther particle reflection until v s h ^ M s c s again. Though 
artificial, this does partly imitate reality, since the shock 
will lose energy as it expands transversely whilst moving 
forward. M s is seen to decrease from 3 to < 2 during 
the simulated shock propagation, suggesting a transition 
between the maximal trapping limit and the isothermal 
limit. By the end of the simulation (t = 95 ps), the ion 
soliton is still visible moving in the plasma with a co- 
herent structure. Therefore, in a cooling plasma v s h con- 
stantly decreases to maintain M = M s , which is itself de- 



creasing. Eventually v s h reduces such that £ m f becomes 
less than the shock thickness, and the shock becomes col- 
lisional. Steepening of the shock profile is observed again 
at very late times due to collisions (figEJ). This phase is 
not modelled in the PIC simulations. 

In conclusion, the dynamics of the shock driven into an 
overdense target due to hole-boring have been observed, 
from the double-layer to ion acoustic soliton and finally 
to collisional phases, for the first time experimentally. 
These results hence provide insight into the dynamics of 
shock acceleration in plasmas, which may prove to be an 
interesting route to high energy ion production. 
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